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ABSTRACT 

We present the results of a ~ 3 year campaign to monitor the low luminosity 
active galactic nucleus (LLAGN) NGC 7213 in the radio (4.8 and 8.4 GHz) and X-ray 
bands (2-10 keV). With a reported X-ray Eddington ratio of 7 x 1CT 4 L Edd , NGC 
7213 can be considered to be comparable to a hard state black hole X-ray binary. 
We show that a weak correlation exists between the X-ray and radio light curves. 
We use the cross-correlation function to calculate a global time lag between events 
in the X-ray and radio bands to be 24 ± 12 days lag (8.4 GHz radio lagging X-ray), 
and 40 ± 13 days lag (4.8 GHz radio lagging X-ray). The radio-radio light curves are 
extremely well correlated with a lag of 20.5 ± 12.9 days (4.8 GHz lagging 8.4 GHz). 
We explore the previously established scaling relationship between core radio and X- 
ray luminosities and black hole mass Lr oc M 6 ~°- 8 L ( ^ 6 , known as the 'fundamental 
plane of black hole activity', and show that NGC 7213 lies very close to the best-fit 
'global' correlation for the plane as one of the most luminous LLAGN. With a large 
number of quasi-simultaneous radio and X-ray observations, we explore for the first 
time the variations of a single AGN with respect to the fundamental plane. Although 
the average radio and X-ray luminosities for NGC 7213 are in good agreement with 
the plane, we show that there is intrinsic scatter with respect to the plane for the 
individual data points. 
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1 INTRODUCTION 

The observable properties of active galactic nuclei (AGN) 
and black hole X-ray binaries (BHXRBs) are consequences 
of accretion on to a black hole at a variety of rates, in a vari- 
ety of 'states', and within a variety of environments. The ma- 
jor difference between the aforementioned classes of object 
is the black hole mass. BHXRBs typically have a black hole 
mass ~1OM while for AGN it is 1O 5 M < M < 1O 1O M . 
Theoretically, the central accretion processes should be rel- 
atively straightforward to scale with mass, and this is sup- 
ported by several observed correlations. These include a re- 
lation between the X-ray and radio luminosities and the 
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black hole mass (Merloni, Heinz & di Matteo 2003; Fal- 
cke, Kording & Markoff 2004), and between X-ray variabil- 
ity timescales, mass accretion rate and mass (McHardy et 
al. 2006). More quantitative similarities between accretion 
'states' and radio jet production have also been demon- 
strated (Kording, Jester & Fender 2006; for the current pic- 
ture of accretion states in BHXRBs and their relation to 
radio jets see Fender, Belloni & Gallo 2004). 

Studying the delays between different emission regions 
gives us a further handle on the scalability of black hole 
accretion, as signals propagate from, for example, the ac- 
cretion flow to the jet. Variability studies have so far shown 
that a correlation exists between the X-ray and optical emit- 
ting regions of both BHXRBs and AGN, typically reporting 
small lags, which are consistent with at least some of the 
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optical variations being due to X-ray heating of the disc 
dRussell et al.ll2009l; iBreedt et alJ feOQQ). A recent study by 



ICasella et al.l ( 2010l ) has shown that a correlated time lag 
of ~ 100 ms exists between the X-ray and IR regions (IR 
lagging X-rays) for the BHXRB GX339-4, indicating a close 
coupling between the hot accretion flow and inner regions 
of the jet. In the case of the BHXRB GRS 1915+105 a vari- 
able X-ray to radio l ag of ~ 30 mins (radio l agging X-ray) 
has been measured l|Poolev fc Fenden 1 19971 : iFender et al.l 
1999). Discrete ejection events have been resolved in both 
the AGN 3C120 iMarscher et alj|2002h and GRS 1915+105 
l|Poolev fc Fenderlll997l : IFender et alJll999h . 

The linear scaling with mass of the characteristic 
timescale around a black hole means that there are advan- 
tages to studying each class of object. In BHXRBs we can 
track complete outburst cycles, from the onset of disc in- 
stabilities through major ejection events, radio-quiet disc- 
dominated states, and a return to quiescence, on humanly- 
observable timescales (typically years). For a typical AGN 
the equivalent cycle may take many millions of years. How- 
ever, for an AGN we are able to resolve individual varia- 
tions of the source on time-scales that are comparable to or 
shorter than the shortest physical time-scales in the system 
(e.g. the dynamical time-scale), something which is currently 
impossible for BHXRBs. In 'black hole time' we are able to 
observe the evolution of sources in fast-forward for BHXRBs 
and in detailed slow-motion for AGN. 

In this paper we present the results of a long term (~ 3 
years) regular monitoring campaign in the X-ray and radio 
bands of the low luminosity active galactic nucleus (LLAGN) 
NGC 7213. Previous X-ray studies show t hat NGC 7213 
is acc reting at a low rate ~ 7 x 10 -4 Lsdd \ Starling et al.l 
2005). The hard state in BHXRBs is typically observed at 
bolometric luminosities below ~ 1% Eddington, and seems 
to be ubiquitously associated with a quasi-steady jet. Above 
~ 1%, sourc es can switch to a softer X-ray state , the jets are 
supp ressed l|Maccarone. Gallo. fc Fenderl 120031 ; iDunn et al.l 
l2009h : furthermore transition to this softer state is usually 
associated with major transient ejection events. As NGC 
7213 is considerably below hEdd ~ 1% we therefore consider 
it a good candidate for comparison with other BHXRBs in 
the low/hard state. If we consider AGN to be 'scaled up' 
versions of BHXRBs by exploring the time lag between the 
X-ray and radio emitting regions we can compare, contrast 
and hopefully relate the accretion and jet production sce- 
narios for AGN and BHXRBs. 

A correlation has been established by ICorbel et al.l 

(2003) and Gallo et al. (2003, 2006) relating the ra- 
dio luminosity (Lr) and X-ray luminosity (Lx) for 
BHXRBs i n the low/hard and quiesc ent states, where 



L R oc L° x 6 . iMerloni. Heinz, fc di Matted d2003h - hereafter 



MHdM03 and iFalcke. Kording. fc MarkofJ (120041 ) extended 
the BHXRB relationship using two samples of AGN to form 
the 'fundamental plane of black hole activity'. By accounting 
for the black hole mass (M) the relationship Lr oc L° x 6 has 
been extended to cover many orders of magnitude in black 
hole mass and luminosity. Further refinements were made to 
the fundamental plane bv lKording. Falcke. fc Corbel! (|2006) 
- hereafter KFC06, using an augmented and updated sample 
to examine the fitting parameters. 

Throughout this paper we define the 'intrinsic' be- 
haviour of AGN and BHXRBs as multiple measurements 



(in the radio and X-ray) of the same source. We define 
the 'global' behaviour as single (or average) measurements 
of multiple sources, both with respect to the fundamental 
plane. 

For the BHXRBs in the low/hard state the relationship 
described above has not only been established globally but 
in some cases intrinsically, i.e GX 339-4, V404 Cyg and a 
small number of other systems have been shown to move 
up and down the co rrelation seen in the fundamental plane 
l|Fender et al.ll2007h . However, in recent years an increasing 
number of outliers have been fou nd below the corr el ation, 
i.e. less radio-loud then expected (|Xue fc Cuil (|2007h : iGallol 
l|2007h ; ICalvelo et al l (|201CD ) as well as some sources which 
move in the plane with a different slope (e.g Ijonker et al.l 
(2010)). To date the correlation found from the fundamen- 
tal plane has only been measured globally for AGN, not 
intrinsically. Note, with respect to the global measurements 
of the AGN population, the specific measurements of the 
radio and X-ray flux used in the correlation are sometimes 
taken at different times and thus could be a source of error 
in the correlation (|Kording. Falcke. fc Corbel|[200r3 ). 

As well as establishing the time lags, another goal of 
this work, is to establish, through quasi-simultaneous ob- 
servations the intrinsic relationship between Lr, Lx and 
M observed in an LLAGN and its relevance to the funda- 
mental plane of black hole activity. We use the MHdM03 
and KFC06 samples for comparison both with an updated 
BHXRB sample taken from Fender, Gallo & Russell (2010) 
- hereafter FGR10. We also explore the possible scatter in 
AGN data points away from the fundamental plane and 
place limits on this deviation. 



2 NGC 7213 BACKGROUND 

NGC 7213 is a face-on Sa galaxy hosting a Seyfert 
1.5 nucleus located at a distance of 25 Mpc assuming 
Ho=71kms -1 Mpc -1 . Narrow low ionisation emission lines 
are observable in its nucl ear spectrum, also m aking it a mem- 
ber of the LINER class l|Hameed et al.ll200ll ). 

The radio properties of NGC 7213 are inter- 
mediate between those of radio-loud and radio-quiet 
AGN. Previous radio studies at 8.4 GHz have not re- 
solved any jet emi ssion from the nucleu s on scales 3 
1 arcsec dBransford et all (|l998T ) , iThean et al.l 



mas to 

(|200ll)lBlank. Harnett, fc Jones! (|2005r i).The Long Baseline 
Array (LBA) is an Australian six station VLBI instrument 
[|Tzioumis![l997h ; LBA observations at 1.4 and 0.843 GHz 
have shown some evidence for l arger scale emission that 
could be due to jet-fed radio lobes ([Blank. Harnett, fc Jones! 
2005). However, as the nucleus remains unresolved it could 
be that the jet is oriented to som e degree in the direction 
of the observer. As suggested by iBlank. Harnett, fc Jones] 
(2005), the lower frequency emission could possibly be 
associated with a 'kink' in the jet at a larger dis- 
tance. T his could be consistent with a gener al model pro- 
posed bv lFalcke. Patnaik. fc Sherwood! (|l996h where radio- 
intermediate objects are in fact radio -quiet objects whose 
jets are to some extent aligned in the observers direction 
and relativistically boosted. 

Previous X-ray studies have so far failed to show signif- 
icant evidence for a soft X-ray excess; Compton reflection; 
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or a broad Fe Kg li n e in the NGC 7 213 sp e ctrum e.g. 
IBianchi et all (|2003h . iBianchi et all <|2008h . IStarling et al 



2005). A narrow Fe Ka line is observed and IBianchi et a" 



2008) show e vidence that it is pro duced in the broad line 



region fBLR). Starling et alj (|2005l ) also report that the ex- 
pected UV bump is either absent or extremely weak. The 
weakness or absence of these signatures suggests that the 
inner accretion disc is missing, perhaps replaced by an ad- 
vection domin a ted a ccretion flow (ADAF) (as suggested by 
IStarling et all l|2005h ) or similar hot flow, consistent with 
the low/hard state interpreta ti on of NGC 7213. 

iBlank. Harnett, fe Jonesl i|2005l ) estimate the black 
hole mass of NGC 72 13 using the nucl e ar ve locity dis- 
persion calculation of iNelson fe Whittle] |l995) and the 
velocity-dispersion v ersus black-hole-mass relationship of 
iFerrarese fe Merrittl (|2000h to obtain 9.6tH X 1O 7 M . 



3 OBSERVATIONS 
3.1 ATCA data analysis 

Bi-weekly monitoring was obtained with the Australia Tele- 
scope Compact Array (ATCA). The interferometer setup 
was such that 128 channels of 1 MHz bandwidth were 
used to form two continuum channels centred at 4.8 
and 8.4 GHz respectively. The radio observations have 
been reduced using the MIRIA D radio reduction package 
|Sault. Teuben. fe Wright! Il995l ). Flux and bandpass cali- 
bration was achieved using (in most cases) PKS J1939-6342 
(B1934-638) and for the phase calibration PKS J2218-5038 
(B2215-508). 

A variety of fitting techniques were then tested to ex- 
tract the flux density of the source from the observations. 
These included testing point and Gaussian fitting to the 
source in the image and uv plane. As ATCA is an East- West 
array and observations were short in duration (~2 hours) an 
elongated synthesised beam was typically produced. Image- 
plane fitting often failed to converge and/or produced non- 
physical results. Therefore fitting was discarded in the image 
plane for the more reliable uv plane method. The final radio 
light curves and spectral indexes (between the radio bands 
only) are shown in figure [T] 

To complement the snapshot observations a full 12 hour 
integration was performed on 2008 June 29; the purpose of 
which was to explore the polarisation properties of the jet, 
and to detect any faint extended structure. Additional steps 
were taken to ensure adequate polarisation calibration was 
performed. This observation was also used to make a deep 
image of the source; beyond the 0.5" beam no weak extended 
radio emission was observed above a noise level a=78/iJ at 
4.8 GHz. 



3.2 RXTE data analysis 

Daily monitoring was obtained using 1 ksec snapshots ob- 
tained with the Rossi X-ray Timing Explorer (RXTE) Pro- 
portional Counter Array (PC A), allowing a long-term light 
curve to be obtained in the 2-10 keV band. The data were re- 
duced using FTOOLS v6.8, using standard extraction meth- 
ods and acceptance criteria. The background was calculated 
from the most recent background models which corrects for 





-4.8 GHz 




°4 




4000 4200 4400 4600 

Time MJD-50000 



Figure 1. Top Panel: X-ray Flux verus Time (MJD-50000). Cen- 
ter Top: Radio flux S v g.4GHz versus time. Center Bottom: Radio 
flux S v 4.sGHz versus time. The errors are too small to show on 
both radio flux figures but are typically ~0.4 mJy/Beam. Bottom: 
Spectral index (8.4-4.8 GHz) a versus time where S 



the recent problems with the RXTE SAA (South Atlantic 
Anomaly) history file. The final 2-10 keV fluxes were calcu- 
lated by fitting a power law to the observed spectra. This 
allows us to take into account changes in the instrumen- 
tal gain over the duration of the monitoring. The final un- 
binned X-ray light curve is shown in figure [T] and is used in 
all subsequent analysis. 



4 RESULTS 

4.1 X-ray/radio light curves and cross-correlation 

Figure [1] shows X-ray flux (top panel), radio flux (middle 
panels) and radio spectral index a where S v oc v +a (bot- 
tom panel) versus time . The X-ray light curve shows a gen- 
eral decrease in flux until MJD 54600. Two distinct X-ray 
flares are observed at MJD 53920 and MJD 54390, both 
appear to be correlated with events in the radio bands. 

We u se the Discrete Correl ation Function (DCF) 
method of lEdelson fe Krolikl | 1988) to calculate the cross- 
correlation coefficients between the entire X-ray and radio 
bands to find the lag. To calculate the centroid lag r cent , 
we use a weighted mean of the positive lags above an 85% 
threshold of the peak DCF value. We use T C ent, rather than 
the peak value r pea k , because the centroid has been shown t o 
better represent the physical lag l|Koratkar fe Gaskelllll99lh . 
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It has been shown that it is difficult to interpret T pea k as a 
physical quantity, and if it is used to calculate the lag, it 
usually offers an und erestimate when compared with T cm t 
|Peterson et alj[l998h . The centroid width is calculated as 



4.8 GHz 



— _L \ " CiTi 
nt ~ N 2-" r. 



(1) 



Where c; and r% are the DCF coefficients and lags, r cen t is 
the centroid lag. As the DCF coefficients are only defined 
for a given lag, we interpolate to find the peak correlation 
coefficient DCF mal at the centroid lag . 

Figure [2] show the DCFs for the X-ray-8.4 GHz 
and X-ray-4.8 GHz bands respectively. We find a lag of 
Tcent =24±12 days at DCF ma:E =0.56 for X-ray-8.4 GHz 
and r cent =40±13 days at DCF maa; =0.7 for X-ray-4.8 GHz. 
The time lags taken from figure [5] are summarised in ta- 
ble [1] The errors in the lag are calculated using the flux 
ra ndomisation / random subset selection (FR/RSS) method 
of iPeterson et afl l|l998h . For 10 4 samples, we randomly re- 
duce the number of data points in our original X-ray and 
radio light curves by 37% and then re-calculate the DCF. 
We take the centroid lag from each of the 10 4 simulations 
and use the rms spread in the distribution of centroid lags 
to obtain the quoted errors. Note that since the radio light 
curves are strongly correlated with one another, the statis- 
tical errors on the lags are not independent between the two 
bands, so that the difference in lag between the two radio 
bands with respect to X-rays is likely to be real (this is also 
highlighted by the radio-radio correlation, see below). 

To ascertain a con fidence level in the DCF c alculation 
we used the method of iTimmer fc Koeniel j 19951 ) to gener- 
ate Monte-Carlo simulations of un-correlated red-noise light 
curves. We assumed a broken power-law-shape power spec- 
trum with low- frequency slope -1, high-frequency slope of 
-2.3 and break frequency v = 0.012 d -1 , based on the best- 
fitting model to fit the X-ray power spectrum (Summons, 
private communication). We used this model to generate 10 4 
random X-ray light curves and cross-correlated them with 
the real radio light curves in both bands. We then used these 
correlations to estimate the distribution of cross-correlation 
values at each lag to obtain a local confidence level, corre- 
sponding to the CCF value which exceeds P per-cent of the 
simulated CCF values at that lag. From the simulations we 
plot on figure 2 the 95% and 99% local confidence levels: we 
therefore assign a local confidence of > 95% for both of the 
peaks with the 4.8 GHz peak reaching close to 99%. 

It is important to note however, that without a priori 
expectation of what the lag should be, it is more statisti- 
cally rigorous to estimate the significance of the correlation 
using the 'global significance', which is the fraction of simu- 
lated random light curves which show CCFs at better than 
the observed confidence level for any given lag. For example, 
to estimate the global significance of a peak which appears 
at the 99% local confidence level, we search the 10 4 simu- 
lated CCFs for any peak values (at any lag) which exceed 
the 99% local confidence level for that lag. The fraction of 
CCFs which do not show such a peak defines the global sig- 
nificance level for the correlation observed at the given local 
confidence level. In this way, we account for the fact that we 
are effectively searching over a broad range in lags, and so 
are sampling from a larger population of potential 'false pos- 
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Figure 2. Top Panel: Discrete cross-correlation function (DCF) 
plot showing time lag against correlation coefficient for X-ray to 
radio (4.8 GHz lagging X-ray) (black crosses). A lag of 40 ± 13 
days at DCF ma:E = 0.7 is calculated using the weighted mean of 
the lags above an 85% threshold of the peak DCF value. Bot- 
tom Panel: Discrete cross-correlation function (DCF) plot show- 
ing time lag against correlation coefficient for X-ray to radio (8.4 
GHz lagging X-ray) (black crosses). A lag of 24 ± 12 days at 
DCFram = 0.56 is calculated using the weighted mean of the 
lags above an 85% threshold of the peak DCF value. On both 
plots the 99% and 95% local significance confidence levels are 
plotted. The errors are calculated via the method described in 
section 4.1 



itives' than we would expect from searching for correlations 
at a single lag. 

We find for the entire range of computed lags (-300 to 
+300 days) the global confidence levels at 4.8 GHz are 82.0% 
(at 99% local significance) and 56.2% (at 95% local signif- 
icance); at 8.4 GHz, 77.1% (at 99% local significance) and 
46.0% (at 95% local significance). If we restrict the range of 
lags to positive only (> days) we find at 4.8 GHz, 91.0% 
(at 99% local significance) and 74.7% (at 95% local signif- 
icance); at 8.4 GHz, 86.6% (at 99% local significance) and 
63.3% (at 95% local significance). 
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Table 1. Summary of the lag times which have been calculated using the Discrete Correlation Function for the X-ray-8.4 GHz and 
X-ray-4.8 GHz light curves. We used the z-transformed discrete correlation function to calculate the lag for the 8.4 GHz-4.8 GHz and 
flare 1 only light curves. 



Time (MJD) 


X-ray-8.4 GHz Lag (days) 


X-ray-4.8 GHz Lag (days) 


8.4 GHz-4.8 GHz Lag (days) 


All 

Flare f (53800-54000) 


24 ± 12 (DCF mM =0.56) 
14 ±11 (DCF maa; =0.78) 


40 ± 13 (DCF ma;c =0.7) 
35 ± 16 (DCF ma:E =0.73) 


20.5 ± 12.9 (ZDCF mQa; =0.81) 



O 0.2 




x 10 



Lag (Days) 

Figure 3. The z-transformed discrete cross-correlation (ZDCF) 
function between the 8.4 GHz and 4.8 GHz light curves(4.8 GHz 
lagging 8.4 GHz). A ccntroid peak of 20.5 ± 12.9 days is found 
at ZDCF m „ x =0.81. 



We note that in the case where we do not restrict the 
lag to lie within any given range, the global significance is 
low. For example, it is not surprising that the 8.4 GHz ver- 
sus X-ray CCF shows two peaks at greater than 95% local 
confidence, since the corresponding global confidence at this 
level is only 56%, so we expect a spurious peak in nearly half 
of the CCFs sampled from random light curves. However the 
peak corresponding to a negative lag of ~ 140 days (i.e. radio 
leading X-ray) cannot easily be explained with any physical 
model. When we presume that there must be a positive lag 
(radio lagging X-ray) - as we would indeed predict based 
on our current understanding of the accretion disk/jet - this 
global significance level increases. Put another way, it is sta- 
tistically unlikely that a spurious correlation should appear 
in the range of lags which matches our physical expectations, 
as seems to be the case here. Therefore, we assume for the 
remainder of the paper that the correlation is real. Note that 
if the correlation is real, the lags can be well-determined: the 
significance of a lag and the significance of a correlation it- 
self are unrelated quantities, since a lag can be very well 
determined from just a single well-sampled 'event' (flare or 
dip) in red noise data, if the data are indeed correlated. We 
will explore possible complexities in the disk/jet connection 
with respect to our calculated lags in section 5.1 



7 - 

\ 6 

'w 5 

at 

E? 

3. 4 
x 

Z3 

E 3 - 

03 

°f 2 
X 

1 






Einstein 1 


A 


Exosat 2 


O 


Ginga 3 


> 


ASCA 4 


< 


XMM-Newton 5 




RXTE 








o' — 

-8000 



1973 
500 r= 



-4000 -2000 2000 

Time (MJD)-50000 



1979 



1984 



1995 



2001 



2006 



2012 



100 



+ 




-8000 -6000 -4000 -2000 2000 4000 6000 

Time (MJD)-50000 

Figure 4. Top Panel: Light curve of archival X-ray 
data points in the 2-10 keV range. Flux values cited in 
ol=|Halpern fc Filippenkd Il984t) A 2 = iTurner fc Pounds! 
119891) , o3 = llNandra fc Pounds! jl994h >4=ASC A using the 
tarturus database ([Turner et al.l l200lh . <5 = llBianchi et al.l 
2008). Bottom Panel: Light curve of archival radio data point s 
at 8.4 GHz. Flux va lues cit ed in Al = dSlee et al.l Il994h . 
< >2=llThean et al.ll200llh o3 = llBransford et al.lll998h and t>4 = 
llBlank. Harnett, fc Jonedl2005h 



4.2 Radio/radio cross-correlation 

To calculate the cross-correlation between the two radio 
bands we used t he z-transf o rmed discrete correlation func- 
tion (ZDCF) of lAlexanderl (|l997n . see figure Note that 
the two radio light curves are sam pled at the same time: as 
noted in lEdelson fc Krolikl l|l988l ): to keep the normalisation 
correct for the DCF we need to omit the zero lag pairs as 
the errors between the two bands are correlated. Therefore 
for convenience (and as the two methods are comparable) 
we switch to the ZDCF which omits zero lagged pairs. We 
calculate the errors using the same method as described in 
section 4.1. We thus find r C ent = 20.5 ± 12.9 days which is 
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consistent with the difference in lags between X-ray-4.8 GHz 
and X-ray-8.4 GHz. 

4.3 X-ray/radio cross-correlation when flaring 

We split the light curve up into data points surrounding 
the first flare at MJD 53920, which is summarised in table 
[1] We then calculated the ZDCF to ascertain an accurate 
measurement of the lag for this flare only. We find T ce nt = 
f4 ± f 1 days for X-ray-8.4 GHz, and T cent — 35 ± 16 days 
for X-ray-4.8 GHz. Both lag times are below the global lag 
for the entire light curve. We do not perform this analysis 
for the second flare at MJD 54380 as it was only sparsely 
sampled in the radio. 

4.4 Long term variability 

The long term X-ray variability, covering a time-range of 
about 30 years, is plotted in the upper panel of figure |U 
Archival data were s earched for in the literature f rom the in- 
struments Einstein (lHalpern fc Filippenkol 19841). EXOSAT 
ijTurner fc Pounds! [j989h . Ginaa (iNandra fc Pounds! 1 19941) . 
ASCA using the T artarus database ([Turner et alj|200l[ ) and 
XMM - Newton dBianchi et al.ll2008l ). Around MJD 44400 
(~ 1980) a sharp flare which is brighter than the flares visi- 
ble in our dataset was observed. Since then the flux has been 
gradually decreasing. 

The long term radio variability at 8.4 GHz is plotted 
in the lower panel of figure [4] The flux values were taken 
from a variety of instruments and publications (which are 
referenced on the plot). In order of in creasing time the 
first point on the plot was taken from ISlee et all (| 19941 ) 
using the Parkes-Tidbindilla interferometer (PTI). Then 
near simultaneous data points wer e taken with the VLA 
|Thean et al.ll200lh and then ATCA jBransford et al.lll998h . 
The flux from the long baseline array (LBA) observation by 
iBlank. Harnett, fc Jonesl l|2005h was the last point plotted 
prior to the ATCA monitoring which is presented in this 
paper. 

When comparing the flux expected from the start 
of our ATCA monitoring with the LBA data point, the 
LBA point seems too low. The core is un-resolved with 
all of the radio telescopes described above. The PTI ob- 
servations constitutes the first data point in figure [4] and 
the LBA the last (before ATCA monitoring). Note, the 
Parkes-Tidbindilla basel i ne is part of the LBA network. 
IBlank. Harnett, fc Jonesl l|2005h comment that in the LBA 
observations, no decrease in flux was observed with respect 
to baseline length i.e it was not resolved; they conclude from 
this that the decrease in flux with respect to the other ob- 
servations is accurate. 

However, both the PTI and the LBA are not sensitive 
to the same spatial resolutions as the VLA and ATCA (the 
largest ATCA baseline is 6km and the shortest LBA baseline 
used in the archival observation is 113km). It is possible that 
some flux (on ATCA equivalent baselines) is not accounted 
for: which could plausibly give a decrease in flux i.e indicat- 
ing brighter, larger scale structure. The logical consequence 
of this is that the PTI data point is also missing some flux as 
it samples on a singular baseline of 275km. We will discuss 
these measurements within the context of the fundamental 
plane and attempt to draw some conclusions in section 5.2 



4.5 Linear polarisation 

A full 12 hour synthesis ATCA observation was per- 
formed on 2008 June 29 at 8.4 and 4.8 GHz. The pur- 
pose of this observation was to perform a reliable polari- 
sation calibration to ascertain an accurate measurement of 
the percentage linear po l arised (LP) flux from the source. 
iBower. Falcke. fc Mellon! l|2002l ) showe d that from a sample 
of 11 LLAGN the mean LP was - 0.2%. lBignall et al.l(l2004l ) 
from a survey of 22 blazars report a mean LP ~ 3% with a 
standard deviation a = 1.5. At the time of observation the 
polarisation of NGC 7213 was < 0.1% at 4.8 and 8.4 GHz re- 
spectively, which is more consistent with the reported value 
of typical LLAGN and not blazars. The percentage polar- 
isation was calculated at other epochs, however without a 
full 12 hour synthesis parallactic angle calibration the values 
were often suspect and will not be presented here. 



5 DISCUSSION 

5.1 X-ray/radio jet connection 

We have shown that a weak but statistically significant delay 
exists between the X-ray and radio emitting regions, with 
the radio lagging behind the X-ray. A number of models have 
been proposed to explain and interpret the X-ray to radio lag 
in both BHXRBs and AG N. The most notable of these are 
the ' internal shock model' dBlandford fc Konigl 19791 ; Ifleesl 
Il978h and a 'plasmon model' ( van der Laanlll966l ). We will 
briefly explore these models and comment on the relevance 
- if any - to our data. 

We first consider a plasmon model. After the accretion 
mechanism(s) have pushed matter into the jet, an adiabati- 
cally expanding initially self-absorbed synchrotron emitting 
plasmon travels at relativistic speeds down the path of the 
jet. At a given time and frequency the matter becomes opti- 
cally thin and is 'detectable' at that given frequency (in this 
case either 8.4 and 4.8 GHz). If the jet is fed at a constant 
mass rate, density and velocity the delay time for material 
to become optically thin will be constant. If these param- 
eters are not constant the delay between X -ray and radio 
will be variable (e.g. see Ivan der Laar] (Il966h) assuming tha t 
the disk and the jet are indeed coupled (|Falcke et al.ll2009h . 
Note, the higher frequency radio emission (8.4 GHz) will be- 
come optically thin first, thus a delay from 8.4 to 4.8 GHz 
is always expected. 

Another model for expl aining the emission seen in j ets is 
the ' internal shock model' (|Blandford fc Konigllll979l ; iRees! 
1 19781 ). The synchrotron lifetime of an emitting region is too 
short to adequately explain the scales of jets observered in 
AGN. These emitting regions - commonly referred to as 
'knots' - are often displaced from the central core emis- 
sion. Localised shocks within these regions are needed to ex- 
plain the time-scales of variation observed (|Reesl 19781 ; iFeltenl 
Il967h . Jet shock scenarios have also been used to model 
the common flat spectrum jet ob s ervered in BHXRBs and 
AGN (e.g. see ISpada et all j|200ll ) ; Ijamil. Fender, fc Kaiserl 
|2009l) ). 

Figure [T] shows the evolution of spectral index; indi- 
cating that during a flare event a increases (a > 0), flat- 
tens and steepens (a < 0) shortly afterwards. The initial 
re-energisation given by a shock would push/compress the 
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Table 2. Summary of archival observations found in the literature. The data points between the horizontal lines indicate where the two 
fluxes were used to calculate the X-ray-radio correlation. Two X-ray and two radio points were taken at similar times thus they were 
both used to calculate the correlation, (R) and (X) denotes the time (MJD) of the observation for radio and X-ray respectively. <5t gives 
the time difference between observations (Radio- X-ray). The instrument used to derive the flux is given in brackets (see section 4.4 for 
references) and errors are given only when reported in the literature. 



Time (MJD) 


Radio Flux S v 8.4 GHz (mjy) 


X-ray Flux 2-10 keV (ergs -1 ) 


Stmeas (Radio - X-ray) (days) 


43985 




2.2X10- 11 (Einstein) 




44168 




1.7X10- 11 (Einstein) 




44374 




3.8xl0~ n (Einstein) 




45644 




5. Ox 10" 11 (Exosat) 




47527 (R) - 48065 (X) 


334±33 (PTI) 


3.5X10" 11 ±0.2 x 10" 11 (Ginga) 


-538 


47527 (R) - 48193 (X) 


334±33 (PTI) 


4.2X10" 11 ±0.25 X lO" 11 (Ginga) 


-666 


49913 (R) - 49479 (X) 


181 ±1 (VLA) 


3.0X10- 11 (ASCA) 


434 


50083 (R) - 49479 (X) 


176 ±1 (ATCA) 


3.0X10- 11 (ASCA) 


604 


51604 (R) - 52057 (X) 


57.6 ±1.3 (LBA) 


2.2X10" 11 (XMM-Newton) 


-453 




0.6 Log + 0.78 Log M 0.6 Log L^ + 0.78 Log M 

Figure 5. Left Panel: The fundamental plane of black hole activity. The large cross indicates the average location of NGC 7213. Black 
dots are the MHdM03 sample with the black line indicating the best fit. Crosses indicate the updated BHXRB sample of FGR10; data 
points have been converted into the 2-10 keV energy range. Right Panel: A close up of the NGC 7213 data points with least squares best 
fit indicated with a dashed line. The MHdM03 fit is shown with a solid line. The archival data points are also shown with but are not 
used in the NGC 7213 data best fit. NGC 7213 measurements are taken at 4.8 GHz and the archival are at 8.4 GHz. 



plasma into the optically thick regime; subsequently this 
process moves the material into the optically thin regime. 
The variable time lag seen in the DCF could be interpreted 
as the time taken for newly injected matter to 'catch up' 
with older matter expanding adiabatically in the jet to shock 
and produce a flare. The lag would be dependent on the 
Lorentz factor of the newly inje cted material whic h in turn 
is related to the accretion rate l|Falcke et al.ll2009T ). 

Although both of these models can adequately explain 
the variations observed in our light curves, we cannot com- 
pletely disentangle them. It is possible that both of these 
models are in some part responsible for the behaviour. As 
the jet remains unresolved we cannot identify an area of 



localised shocks or indeed discrete resolved events that we 
could use to support a particular model. 

To explore - in a very simplistic sense - how the time 
lags seen in BHXRBs scale up to AGN we offer two com- 
parisons; that of a simple mass scaling and also mass- 
Eddington ratio scaling. For example, discrete resolved ejec- 
tion events have been seen in GRS 1915+105 with a variable 
time lag between X-ray and radio (GHz), with the clear- 
est ex amples of events having lags of Tgrsiqis = 20-30 
mins (|Poolev fc Fender! Il997l : Mirabel et al.lll99Sl ). Taking 
the characteristic timescale of X-ray-radio lags measured in 
GRS 1915+105 and scaling up to NGC 7213 with mass 
only (assuming the mass of GRS 1915+105 as ~ IOMq 
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25 I 1 1 1 1 1 1 

18 20 22 24 26 28 30 32 

0.63 Log L + 0.64 Log M 



Figure 6. The fundamental plane using the KFC06 LLAGN sam- 
ple only (circles) and updated FGR10 BHXRB sample (crosses). 
The large cross indicates the average location of NGC 7213. The 
solid line indicates the KFC06 best fit. 

and the mass of NGC 7213 as 9.6 xl0 7 M Q and using 
Tiag = Tgrsi915 x ( M N GC72i3 /Mgrsi915 )) we infer T iag = 2 
xlO days, much longer than we observe. GRS 1915+105 is 
however accreting at a much higher rate, therefore including 
a scaling by the ratio of the Eddington luminosities (using 
Tiag2 = Tiag x {L E dd-NGC72iz/ L E dd-GRSWi5) and taking 
L E dd-GRSi9i5 ~ 1) we find Ti ag2 = 140 days, reasonably 
comparable to our measured lag. 

The difficulty in comparing the actual X-ray-radio time 
lag measured for NGC 7213, and that of the GRS 1915+105 
scaled time lag is that the radio emitting regions being 
probed are significantly different: in a broader view the spec- 
tral energy distributions are different. In fact until the struc- 
ture of jets and how it scales with accretion rate and mass 
at a given frequency are better understood, such attempts 
at quantitative comparison are of marginal value, although 
the qualitative comparison is valuable. 

5.2 The fundamental plane of black hole activity 

The fundam ental plane of black hole ac- 
tivity dMerloni. Heinz, fc di Matted 120031 ; 

iFalcke. Kording. fc Markofj |2004 ) showsthat the cor- 
relation Lr oc M 0S L X 6 holds over many orders of 
magnitude of X-ray and radio luminosities, and black hole 
masses. Figure [S] shows; at high luminosities a sample of 
AGN taken from MHdM03, and at the low luminosities 
an updated BHXRB sample taken from FGR10. The 
updated BHXRB X-ray data were observed in the 0.5-10 
keV energy range, therefore careful steps were taken (using 
WebPIMMS) to correctly convert into the 2-10 keV range 
to compare with the MHdM03 and NGC 7213 data points. 
We plot on figure[5]the best fit parameters Lr = (0.6jlo 11) 
log L x + (0.78±o.o9) log M + 7.33±l° 5 7 as defined by 
MHdM03. Note, we do not re-fit the best fit line for the 
updated BHXRB sample. 

We paired the X-ray and radio data points to calcu- 



late the correlation by finding the closest X-ray point (in 
time), to the radio, because the X-ray sampling was more 
frequent. The average for NGC 7213 sits well on the pre- 
dicted correlation. The right hand panel of figure [5] shows a 
close up of the NGC 7213 data. A least squared best fit for 
the NGC 7213 data points only is shown and parametrised 
by Lr = (0.2±g:g||) log L x + 28.3±\'P 7 (the constant in- 
cludes the mass term in the context of MHdM03). 

We include in figure [5] the correlation found from the 
archival X-ray and radio data; however, as they were not si- 
multaneous we do not include them in the best fit. We paired 
the X-ray and radio data points to calculate the correlation 
according to the closest in time. If either two X-ray, or two 
radio points were very close in time, both are included in 
the correlation. Also note that these data were calculated 
between 8.4 GHz and 2-10 keV and not 4.8 GHz. Assuming 
that the same correlation (i.e ~ flat spectrum) holds be- 
tween the 8.4 and 4.8 GHz observations the plot seems to 
show that two of the archival points do not follow the trend 
predicted by the fundamental plane. 

To assess whether the deviation/scatter in the archival 
points (and indeed the other MHdM03 AGN) away from the 
fundamental plane can be explained due to a delay in sam- 
pling the X-ray and radio fluxes; we take the fraction RMS 
scatter in the RXTE/ATCA X-ray and radio data and plot 
this errorbar with the archival points. We find the fractional 
X-ray RMS to be fRMS x =26% and the fractional Radio 
RMS fRMS r =ll%. The full light curve presented in this 
paper spans ~ 1000 days and the longest change in time 
between the archival observations is ~ 500 days (see table 
[5}. Although this method only offers an estimate of the er- 
ror between sampled data points, it does allow us to assess 
outliers. 

Within this framework it appears that the first (PTI 
correlation) and last (LBA correlation) data points that 
were used in the correlation are outliers: the ATCA and 
VLA correlation points sit close to the best fit. Even when 
taking into account the full RMS from the 1000 day light 
curve the error bars do not bring the data points close to 
the MHdM03 best fit. We have also speculated earlier that 
there could be a certain amount of missing flux associated 
with these long baseline observations. Moving both the PTI 
and LBA points up by some set amount will still leave the 
PTI point away from the correlation. Although applying the 
intrinsic variation of the source to the archival points cannot 
bring all points onto or close to the correlation; the scatter in 
these points is within that permitted by the other MHdM03 
points shown on the plot. 

By using the fractional RMS of X-ray and radio vari- 
ability from our study we have placed a constraint in the 
deviation of the archival data points away from the funda- 
mental plane. These constraints suggest that other forms of 
scatter (apart from bad sampling and missing flux) could 
be affecting the data. As summarised in KFC06 other forms 
of scatter could be attributed to beaming, source peculiari- 
ties and spectral energy distribution - but note, these should 
remain relatively constant for the same source. 

In figure [6] we plot a LLAGN only sample and best fit 
parameters taken from KFC06 with the updated BHXRB 
sample. Note, the X-ray data in the KFC06 sample are taken 
in the 1-10 keV band. As the NGC 7213 data points are in 
the 2-10 keV band we assume a photon index of r=1.8 and 
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add a correction factor to the NGC 7213 data points to 
make them comparable. We apply a similar correction to 
the BHXRB sample. From this plot it is clear that the NGC 
7213 data points are positioned slightly above the best fit 
line. 

Considering the postion of NGC 7213 on the funda- 
mental plane with respect to other LLAGN (see figure 6), 
we calculate its radio loudness parameter to assess the 
differentiation. We calculate the radio loudness parameter 
R = L(j cm / Lb (where La cm and Lb are the radio and opti- 
cal luminosities); we use a B band magnitude of 16.3, and 
find the optical flux So V t u sing B — —2.5\ogS op t — 48.6 
l|Halpern fc Filippenkol 1 19841 ): giving R = 134.8. In this 
scheme radio-loud sources are typically defined as having 
an R parameter > 10, w hile radio-quite range between 0.1-1 
l|Kellermann et al.|[l989h . 

Using the alter n ative radio loudness parameter of 
iTerashima fc Wilson! (2003) which utilises the X-ray in- 
stead of optical luminosity, R x = L ecm / L 2 - iokeV , we find 
logRx = -3.28. iPanessa et al.l l|2007h show that for a sam- 
ple of low-luminosity Seyfert Galaxies logi?x = — 3.64±0.16 
while for a sample of low-luminosity radio Galaxies (LLRGs) 
logRx = —1.40 ±0.11. Therefore, with respect to the X-ray 
radio loudness, NGC 7213 is only slightly higher than that 
of a sample of low-luminosity Seyfert Galaxies; while under 
the standard definition of radio loudness it is indeed radio 
loud. These results are consistent with the position of NGC 
7213 on figure H 

As was discussed earlier in this paper there is an appar- 
ent time lag between events in the X-ray and radio. There- 
fore comparing the fitting parameters found from the NGC 
7213 data with the MHdM03 relationship without correcting 
for the lag might give rise to errors as we are not matching 
the correct data points. The width in the cross-correlation 
peaks shows that the time lag is variable. Thus, for example, 
the two radio flares could have different lag times associated 
with them. Therefore shifting the entire radio light curve 
back by a set amount to match the X-ray could still give a 
scatter. To simplify this problem we separated out the data 
for the first flare only because we have a more accurate mea- 
surement of the lag in this specific case. We then shifted the 
radio data -35 days which was the time lag measured for 
this singular flare using the DCF at 4.8 GHz (see table [T]). 
The top panel of figure [7] shows the uncorrected data on the 
MHdM03 plot while the middle panel shows the corrected 
data; for completeness the bottom panel shows all radio data 
points shifted back. 

For the first flare correcting the data appears to reduce 
the scatter and increase the gradient more in line with the 
MHdM03 best fit. It is now described by L R = (0.58t°:w) 
log Lx +8.5^£'g. To check the statistical significance of this 
we measured the gradient for a variety of shifts. From 0-25 
days the gradient gradually steepens until it gets close to 1 
(giving a coefficient of ~ 0.6). From 30-50 days it plateaus ~ 
1 and from 50 days the gradient decreases towards 0. Thus 
it does appear that moving the flare back by the amount 
given from the DCF function does seem to better represent 
the data with respect to the MHdM03 fit. However, it should 
be noted that as this is a log/log plot, measuring gradients 
from such a small range of values should be treated with 
care. It would, however, be interesting in future studies to 
assess the importance of this parameter. 
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Figure 7. Top Panel: Flare 1 radio data points at 4.8 GHz un- 
corrected for lag. On both graphs the solid line indicates the least 
squared best fit for the NGC 7213 data only; the dashed line in- 
dicates the MHdM03 best fit. Middle Panel: Flare 1 radio data 
points shifted -35 days to correct for time lag calculated from the 
cross-correlation function. Bottom Panel: All radio data points 
shifted back -35 days. 



6 CONCLUSION 

We have used the Australian Telescope Compact Array and 
the Rossi X-ray Timing Explorer to conduct a long term 
study of AGN variability in the LLAGN NGC 7213. We have 
used the cross-correlation function to show that a complex 
and weakly significant correlated behaviour exists between 
the X-ray and radio emitting regions. Although the statistics 
only show a weak correlation, this study is the first definitive 
campaign to probe this type of behaviour. 

We have shown that NGC 7213 sits well on the pre- 
dicted fundamental plane of black hole activity plot when 
compared with the MHdM03 sample. However we have 
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shown that when comparing NGC 7213 with a revised 
BHXRB and LLAGN sample that the data points are above 
the expected correlation; which is however consistent with 
the calculated radio and X-ray loudness parameters. We 
have also shown some support that by correcting for the 
time lag between events in X-ray and radio the gradient of 
data points agree better with the best fit derived from the 
MHdM03 sample on the fundamental plane. 
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